Period-doubling-bifurcation readout for a Josephson qubit 
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We propose a threshold detector with an operation principle, based on a parametric period- 
doubling bifurcation in an externally pumped nonlinear resonance circuit. The ac-driven resonance 
circuit includes a dc-current-biased Josephson junction ensuring parametric frequency conversion 
(period-doubling bifurcation) due to its quadratic nonlinearity. A sharp onset of oscillations at 
the half- frequency of the drive allows for detection of small variations of an effective inductance 
and, therefore, the read-out of the quantum state of a coupled Josephson qubit. The bifurcation 
characteristics of this circuit are compared with those of the conventional Josephson bifurcation 
amplifier, and its possible advantages are discussed. 
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The problem of an efficient readout of solid state quan- 
tum systems including Josephson qubits (see, e.g., Ref.Hl) 
is of high importance from both theoretical and practical 
points of view. The dispersive readout techniques based 
on the radio- frequency measurement of reactive electrical 
parameters (for example, the Josephson inductance^ or 
quantum Bloch capacitance 3 -^) received significant recog- 
nition, since they allow one to minimize the backaction of 
the readout circuit on a Josephson qubit. Recently, par- 
ticular interest was focused on such systems operating 
in the non-linear resonance regime (Duffing oscillator), 
which was possible due to a cubic non-linearity of the su- 
percurrent in a zero-phase biased Josephson junction^— . 
In this regime, under the action of a weak signal and/or 
fluctuations the circuit undergoes a bifurcation, i.e., a 
transition between two stable oscillatory states^. The 
successful idea of application of such a Josephson Bifur- 
cation Amplifier (JBA) for measurements of a qubit was 
first proposed by Siddiqi et alfiS, and this has served 
for us as a motivation for the development of a readout 
based on another type of bifurcation in superconducting 
non-linear circuits. 

In this paper we propose a readout circuit, whose op- 
eration principle is based on excitation of half-harmonic 
oscillations, i.e., a Period Doubling Bifurcation (PDB), 
occurring in an rf-driven resonance circuit with quadratic 
non-linearity of reactance. The essential difference of the 
PDB from the JBA regime consists in the parametric na- 
ture of the PDB resonance manifesting itself in abrupt 
switching from zero-oscillation state into the dynamic 
state with a double period and appreciable amplitude 
of the oscillations 1 ^. This regime may be favorable for 
an output-stage preamplifier receiving in the case of the 
PDB a signal with zero background. Moreover, as we 
shall show below, the switching characteristics of our cir- 
cuits are somewhat different from those of conventional 
JBA; in particular, we find that in addition to better 
contrast between two possible stationary states of the 
PDBA, it may have a narrower switching region. 

The PDB circuit (see Fig. [T]) comprises a dc-current- 



biased Josephson junction with the critical current I c , ca- 
pacitance C including the self-capacitance of the junction 
with, possibly, a contribution of an external capacitance, 
the linear shunting conductance G, as well as an attached 
qubit, presented here as a charge-phase qubit 2 -- 2 . The 
circuit is driven by a harmonic signal J ac = I a cos 2ujt 
at a frequency close to the double frequency of small- 
amplitude plasma oscillations oj p , i.e., u ~ lo p . 

Neglecting fluctuations, the dynamics of the bare sys- 
tem (excluding the qubit, whose quantum state only 
slightly changes the plasma frequency of the entire cir- 
cuit, uj p — > uip) is governed by the model of a resistively 
shunted junctioni 3 -: 
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where the finite current bias Iq < I c ensures a dc phase 
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FIG. 1: Electric circuit diagram of the period-doubling bifur- 
cation detector with microwave-based readout. The resonator 
is formed by the inductance of a non-linear Josephson junc- 
tion (large crossed box), biased at a non-zero phase value tpo, 
and the capacitance C. 

The linear losses are accounted for by the conductance G. The 
resonator is coupled to a charge-phase qubit formed by a su- 
perconducting single electron transistor with capacitive gate 
(left) and attached to the Josephson junction. The qubit op- 
eration at the optimal point for an arbitrary bias Jo is ensured 
by a proper value of the external magnetic control flux $ c , 
applied to the qubit loop, and the gate charge Q g on the qubit 
island. 
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drop ipo = arcsin(/ /P) across the Josephson junction. 
The small-ac-signal expansion {x -C 1) of the Joseph- 
son supercurrent term includes the following components: 
sintys = sm{<po+x) ~ sin<po(l — x 2 /2+x 4 /24)+cos(^o(^— 
x 3 /6). The angular frequency of small oscillations of ip 
around tpo is uj p = (cos (/?o) 1 ^ 2 w p o, where the bare plasma 
frequency is lu p q — (2e/ c /?iC) 1 / 2 . 

Using the dot to denote derivatives with respect to 
the dimensionless time r = uit, we write the equation of 
motion for x in the form: 

x + x = £x ~ 26x + /3x 2 + 7X 3 - nx 4 + 3P cos 2r. (2) 

The dimensionless coefficients in this equation are 

£ = 1-k<1, 6 = G/2uC = 1/2Q < 1, (3) 

/3 = = (/stan^o)/2, 7 = /c/6, 3P = kIa/I c , 

where ft = (lo p /uj) 2 ~ 1 and Q is the quality factor. 
The quadratic non-linear term (oc /3) in Eq. @ ensures 
parametric down-conversion from the drive frequency 
2cj. Note, that similar frequency conversion, and there- 
fore the PDB effect, is also possible without quadratic 
non-linearity, if instead of the driving force 3P cos 2r in 
Eq. (|2|) the circuit is parametrically driven by a term of 
the form oP'cos2t (see, for example, Refs. [lllflih . This 
case can be realized by a periodic modulation of the crit- 
ical current of the zero dc-biased Josephson element I c in 
a split dc-SQUID configuration by using an alternating 
magnetic flux driving in the SQUID loopi£. 

The leading terms in the solution of Eq. © have 
the form x = y — Pcos2r, where y(r) denotes oscilla- 
tions at the frequency s» u, and the second term is the 
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FIG. 2: (Color online) Intensity A 2 of oscillations of the 
Josephson phase at frequency cj versus frequency detuning 
£ for two amplitudes of the pumping signal (at frequency 
2b)). Dashed lines show unstable states. When the detun- 
ing approaches a bifurcation point ((i) or (ii)), PDB occurs 
(vertical arrows from or to the zero state, respectively). For 
comparison, a typical resonance curve of a JBA is sketched in 
the inset. 



forced oscillation at the drive frequency 2uj; x(t) also 
contains other harmonics at multiple frequencies, which 
strongly influence its dynamics and stationary states^, 
cf. Eqs. (fTTI TT^|) below. We apply the method of slowly- 
varying amplitudes by introducing slow variables^, 



y = A cos(r — a) , y rs —A sin(r — a) . 



(5) 



The variables A(t) and cx(t) present the amplitude and 
phase (relative to the drive) of the oscillation at the half- 
frequency of the drive; they vary weakly over the period 
of these oscillations (with dimensionless rates <C 1). Ac- 
cordingly, 



= A cos ( 



v = A sin ( 



are two quadratures of these oscillations, u 2 



(6) 

v 2 = A 2 . 



The dynamics of the slow variables is governed by the 
equations: 



.4 

Aa 



T+2jr dr . / sin(r - a) 

— eF(x, x, t) ' 



2ir 



cos(t — a) 



(7) 

with averaging over a 2-7r-period of the oscillations at fre- 
quency cj (1 in dimensionless units), where the function 
in the integrand 



eF(x, x, t) =£x- 26x + 3x 2 + 7a: 3 - \ix 4 



(8) 



includes small terms at frequency lo and large terms at 
the drive frequency 2u and its higher harmonics. The 
averaging over the period of oscillations in Eq. ([7]) yields 
a pair of reduced equations for the amplitude and phase: 

A = -6A- -Asm2a(j3P - fiPA 2 ) , (9) 

£ BP 3 
a = - - — cos 2a + -jA 2 + jlPA 2 cos 2a . (10) 



The coefficients f(P), /3(P), j(P), jl(P) to the leading 
order in P 2 are given byi£ 



fj, 



8 = P, 
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7 = 7 + y 



Ufa +1.0*. 
P 16 ' 2<T 



(11) 
(12) 



which implies that 7 « (3 + 5 tan 2 (p )/18 and jl « 
(7/192) tan ip (tan 2 ip — 1). Corrections of order P 2 to 
these coefficients do not change further analysis qualita- 
tively, but only slightly modify the results quantitatively. 

Equation ([§]) always has a trivial solution A = 0. In 
the limit of weak pumping (Pel) and small resulting 
oscillations (A 2 <C 1), the last terms (oc jlP) on the right 
hand side of Eqs. (|9T TTU|) can be neglected, and the os- 
cillation amplitude of the non-zero stationary solutions 
(A = a — 0, A 7^ 0) may be found explicitly-ii: 



A± = — 

± 37 



£±J{/3P) 2 -49 2 



(13) 
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FIG. 3: (Color online) Stability diagram of the period- 
doubling bifurcation readout. Stationary oscillatory solutions 
appear and disappear as the detuning £ and driving ampli- 
tude P are varied. This is indicated by the sketched potential 
curves, which show stable (minima) and unstable (maxima) 
states. The former include: only the zero solution in the lower 
region, only a pair of equal-amplitude solutions A+ , A* + with 
a 7r phase shift in the upper region, and the zero and the 
pair A+, A* + in the left, 'bistable' sector-shaped region. The 
equations in the figure describe the bifurcation boundaries. 



For the pumping amplitude exceeding the threshold set 
by dissipation, \f3P\ > 28, the values 



e± = ±d(pp) 2 -w 2 



(14) 



yield the range of frequency detunings, £_ < £ < £+, 
within which the zero solution is unstable. In this range 
the system switches into the oscillating state with a fi- 
nite amplitude A + given by Eq. (1131) . For £ < £_ the 
parametric resonance curve is multivalued with the sta- 
ble trivial A = and nontrivial A + solutions, while the 
solution A_ is unstable. Taking into account higher (e.g., 
oc p) terms in Eqs. (I MTU)) ensures that -<4+(£) and A-{£) 
merge, limiting both the amplitude A + and the the range 
of bistability in £; for stronger drive even higher nonlin- 
earities become important. The shape of the resonance 
curve, calculated numerically from Eqs. (|9lHU|. is shown 
in Fig. [2] for several values of the drive amplitude 3P 
just above the excitation threshold. However, for fur- 
ther considerations of the threshold behavior the higher 
nonlinearities are not crucial, and below we neglect the 
/i-terms. 

The stability diagram of the system in the space of 
the control parameters, the detuning £ and the driving 
amplitude P is shown in Fig. [3] The parameter plane is 
divided in three regions, with the following stable-state 
amplitudes (cf. Eq. ([TB")) ): A = in the lower region, A + 
in the upper region, and both and A + in the 'bistable' 
sector (this region is limited by two solid lines). The 
bifurcation lines are given by the relations A + = A_ 
(lower left horizontal solid line), A_ = (i.e., £ = £_(P), 



upper solid curve) and A + — (i.e., £ = £+(P), dashed 
curve) . The coordinates of the triple point are £ = and 
P = 26/(3 fa 2/(<2tan<p )- 

Equations for the quadrature components of the veloc- 
ity field, 



A cos a — Ad sin a 
A sin a + Ad cos a 



(15) 



where A and d are given by (r31 ll0[) . can be represented 
as Hamiltonian equations of motion with friction: 



-d v H — 9u, 



d u H - 9v , 



or, equivalently, 



(16) 



(17) 



A = -6A - -d a H, Ad = d A H , 
A 

where the Hamiltonian is given by: 



H=(t-pPcos2a)^- + ^A\ (18) 

This Hamiltonian for the slow variables can be obtained 
from the Hamiltonian for the physical quantities^. Fig- 
ure S] shows a contour plot of the absolute value of the 
velocity |v| = (ii 2 + v 2 ) 1 / 2 = {A 2 + A 2 d 2 ) 1/2 in the case 
of a multivalued stationary solution. One can see the 
darker S'-shaped narrow valley, where the motion is slow 
along the curvilinear s-axis. The black spots in this area 
show the stationary solutions, which are the stable focus 
at zero, A — 0, the stable foci A + and A*, correspond- 
ing to equal-amplitude oscillations with a mutual phase 
shift of 7T, and the unstable saddles A- and A*_ (also 
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FIG. 4: Contour velocity plot calculated from Eqs. (f MTO)) 
and (|15[) for the parameters corresponding to a multivalued 
solution. The absolute value of velocity |v| is lower in darker 
areas. The curvilinear trajectory along the valley of minimal 
velocity (short-dash line) is parametrized by s. The straight 
dashed line at origin indicates the most probable direction of 
escape from the zero state close to the bifurcation. 
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with a mutual 7r-shift). For weak dissipation these 'sad- 
dle points' are the lower points of the barriers separating 
the basins of attraction of the foci in the landscape of H. 
Thus the most probable escape path from the zero state 
is along the S'-shaped valley. 

In the vicinity of the bifurcation point £_ within the 
bistable region, the height of the energy barriers is small, 
and one can show that there is a separation of time scales, 
which can be used to solve the dynamics: the fast relax- 
ation from outside towards the .S'-shaped valley is fol- 
lowed by slow dynamics along the valley. In this region 
the points A_ and A*_ are close to the origin, u = v = 0, 
and the slope ao of the valley at the origin can be found 
from Eq. 

sin2a = -?7~\ where r\ = (3P/26. (19) 

To describe the slow motion along the valley near the 
origin, where one can use the amplitude A as a coordi- 
nate, we first solve an equation for the fast motion in the 
axial a-direction (variable a relaxes fast, with a typical 
rate of 9). To find the subleading nonlinear terms in the 
equation of motion along the valley, one needs to take 
into account the deviation of the valley near the origin 
from a straight line. The resulting equation of motion 
can be represented in the form of an easily solvable ID 
equation (cf. Ref. 0): 



ds 



dW(s) 



(20) 



where the pseudopotential is to the lowest orders a bi- 
quadratic polynomial, W = as 2 — 6s 4 , where 



a = - — y/r] 2 - 1, 



b = — W ?i 2 — 1 . 
32 IV 1 



(21) 



For £ < £_ (5£ < 0) one finds that a > 0. Thus, when £ 
crosses £_ from above, the zero unstable stationary solu- 
tion bifurcates and separates into a stable solution at zero 
and two symmetric unstable solutions A_ (Fig. [5]). This 
property makes it sensitive to small changes of the circuit 
parameters (in particular, to the qubit state via its effec- 
tive Josephson inductance, which modifies the detuning 
£). The switching characteristics of such a detector can 
be found from the analysis of this system in the presence 
of noise, which results in a finite width of the transition. 
To describe the bifurcation-based readout, one needs to 
find the tunneling rate out of the shallow well W(s) near 
the bifurcation. 

Small fluctuations due to the conductance G are taken 
into account by adding a noise term SI, with the spectral 



density Si(ui) = 2?kjGcoth 



2/c B T> 



to the rhs of Eq. (JTJ. 



This gives rise to independent fluctuations of the two 
quadratures. Their correlation functions are 

(Su(t)Su(t')) = (Sv(t)Sv(t')) = 2T cS S(r - r'), (22) 

with (5u(t)Sv(t')) = or 

(6A(t)SA(t')) =A 2 {6a( T )8a(T')) = 2T cS 8(t-t'), (23) 
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FIG. 5: A sketch of the potential W(s) along the S'-shaped 
region (see Fig. [4]) for two ultimate cases: (a) bifurcation 
— » A + (AD (near the upper solid line in Fig. [3]) and (b) 
bifurcation A+ — > (near the lower solid line in Fig. [3} . Near 
the bifurcations the dashed lines approximate the energy bar- 
riers by a quartic and cubic polynomials, respectively. 



with (6A(T)Sa(T')) = 0, where the effective temperature 
k 2 Si(uj)uj low-ui K 2 fcsTGw 



T r (f = 



8(/ c cos ip ) 



2(/ c cos c/? ) s 



(24) 



and the latter expression holds in the low-frequency (clas- 
sical) limit Tvjj -C fegT. Upon the reduction to the ID 
equation ([2"fjf this gives fluctuations with the same noise 
power, 



(8s(t)5s(t')) = 2T eS S(i 



(25) 



which affect the motion along the s-coordinate. 

Adding the Langevin term 6s(t) on the right-hand side 
of Eq. (|20[) . one can derive and then solve a ID Fokker- 
Planck equation^ (in fact, a Smoluchowski equation since 
the 'mass' term, oc d 2 s/dT 2 , is absent in Eq. ([20])) for the 
probability density w(s, t) to find the system at the point 
s at the time t, 



dw 

at 



d_ 



8W(s) 

— 1 

OS 



'eff 



ds 2 



(26) 



The escape rate out of the zero metastable state A = 
to the stable state A + or A* + is given by the Kramers for- 
mula— reflecting the activational behavior of the system, 



T = 2(u A /2n)e- Aw ^ 



(27) 
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FIG. 6: (Color online) Switching curves for the period- 
doubling bifurcation readout at various temperatures. Note 
saturation at low T, governed by 'quantum noise' at fcsT < 
huj. Inset: low-T saturation, illustrated by the curves' width 
<5£sw (inverse slope at P = 0.5, oc y/S(u))) as a function of 
6 = y/k B T/hui. 



where the factor 2 accounts for two escape possibilities 
(to the left or right wells). For the overdamped case 
of a zero-mass particle the formula for uja is given, for 
example, in Ref. Il9l . The pref actor uja is determined by 
the geometrical mean of the curvatures of W(s) at the 
bottom of the central well (equal to 2a) and at the top 
of the barrier (equal to 4a), i.e., 



UJA 



2\/2aw = uj^y/r) 2 - 1. 
\/2 



(28) 



The barrier height to the lowest order in d£ = £ — £_ is 
AW = a 2 /4boc (£~£-) 2 - (29) 

Typical switching curves (switching probability P sw — 
1 — e~ rr ° during some observation time t vs. £) are 
shown in Fig. [6] for various temperatures for a set of typ- 
ical circuit parameters. Note that the position and the 
width of the switching curve (see inset) saturates at low 
temperatures. This effect is not a manifestation of the 
real quantum tunneling, but is rather linked to the fact 
that activation in the rotating frame of the first harmonic 
Eqs. (|5I6[) . i.e., the low- frequency noise in that frame, is 
given in the laboratory frame by the noise at a finite 
frequency u, cf. Eq. (|24l) and above. 

Equation ([29)) implies that the width (along the detun- 
ing axis £) of the switching curves, given by the inverse 
slope (dP/d£) _1 at P = 0.5, scales as (5£ sw oc y/T above 
saturation. Thus as T — > 0, it falls off slightly slower 
than that for the 'standard' Josephson bifurcation am- 
plifier—^, where 5£ sw oc T 2 / 3 . This (minor) difference 
stems from symmetry of the PDBA w.r.t. a shift by a 
drive period: a — > a + 7r. This symmetry implies that the 
generic form of the ID potential W near the bifurcation 



is es 2 — as 4 unlike es — as 3 for the JBA. Here e measures 
the distance from the bifurcation, and s is the relevant co- 
ordinate in phase space. However, this symmetry can be 
broken, and the stronger effect of cooling (S^ sw oc T 2 / 3 ) 
restored, by a weak admixture at frequency ui to the drive 
signal. An alternative strategy consists in using another 
bifurcation point, where A+ = A- in Fig. [5] (on lowering 
the detuning £, the system follows the solution A + until 
it merges with A_, where it switches abruptly to zero; 
there is no symmetry around this point). 

Thus we have suggested two protocols of operation of 
the PDBA (with potentials shown in Fig. [5] and oper- 
ation indicated by arrows in the parametric-resonance 
plot Fig. [5]): one of them involves switching from the 
zero state to a large-amplitude stable state A + near the 
bifurcation point and the other involves a reverse 
switching from the large-amplitude state A + to zero near 
the merging point of A + and A_ . Note that in both cases 
to perform a read-out, that is to find out if a switching 
has occurred, one needs to distinguish a zero state from 
a large-amplitude state. This should be contrasted with 
the JBA, where two finite-amplitude (and often, similar- 
amplitude, but different-phase) states have to be distin- 
guished. From this viewpoint, the PDBA may be more 
convenient in practical applications. Other protocols can 
also be discussed (cf. Ref. H^) ■ 

The readout of a coupled qubit is based on the shift in 
the plasma frequency (and thus of the switching curve) 
due to different Josephson inductances in two qubit 
states. The inductance values depend on the type of 
qubit and its parameters. During the readout the qubit 
state is encoded in the resulting oscillations of the PDBA 
by tuning the control parameters (such as the drive fre- 
quency and amplitude, i.e., £ and P) to a point with the 
maximal difference (contrast) between the two switch- 
ing curves. High contrast is reached when the shift in 
the plasma frequency exceeds the width of the switch- 
ing curve. In an ideal arrangement, this contrast reaches 
100%: P sw = and 1 for two qubit states. For the PDBA 
the contrast reaches values comparable to those for the 
JBA with similar circuit parameters (for example, about 
0.3% in frequency sensitivity for the parameters of Fig. [5] 
at low T, that is sufficient for reliable readout of the 
charge-phase qubit — shown in Fig.Q]). Further optimiza- 
tion of the PDBA parameters is possible. 

Let us compare the switching curves for the PDBA and 
JBA 5 near the upper critical lines of the bistability region 
(iff in the notations of Ref. @) . We consider the tunneling 
exponents as functions of the dimensionless deviation of 
the drive amplitude from the bifurcation, 1 — P j P- , for 
the PDBA, and we use the same notation, instead of 
1 — rj/rj-f, for the JBA. According to Refs. HJH, for the 
JBA 



1 



P 

~pZ 



Ej\t\ 



2/3 



(30) 



For the PDBA, we replace the difference (£ — £_) in 
Eq. (|2"9"j) by 9p£_(P_ — P). Assuming that the detuning 
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£ > 1 (that is 5uj > u/Q), we find that the tunneling 
exponent is of order 



AW ^2^7 cos ip E j 



cir 



k B T 



P 



(31) 



Thus increasing the detuning (and the corresponding 
driving amplitude P_(£)) can suppress the width of the 
relevant switching curve (switching probability- vs. -drive 
amplitude P): 



P 

~p~l 



k n T 1 



1 



£/|?|Q 2 |£| 2 6 7 cos^ 



1/2 



(32) 



We note that various operation protocols of the read- 
out device based on PDBA are possible, and one can force 
crossing the bifurcation region and switching between the 
oscillating states by tuning various parameters. In par- 
ticular, the current bias Iq, the amplitude and frequency 
of the drive can be used for engineering a metapotential 
of desired shape and, therefore, optimization of the read- 
out. In our analysis we have focused on the noise-induced 
activation over the barrier in this metapotential. As we 
saw (cf. above Eq. ([2"3|)). the effective 'temperature' is set 
by the noise level at frequency u> and saturates on lower- 
ing the temperature T below huj/ks- This low-T regime 
may also be thought of as 'quantum activation ' 20 ! 21 . One 
could also consider the quantum tunneling^. However, 



in similar systems the corresponding tunneling rate is ex- 
ponentiallysmall, especially close to the bifurcation point 
(cf. Refs.SH. 

In conclusion, we have suggested to use a nonlinear 
Josephson resonator, driven near its double plasma fre- 
quency, as a sensitive quantum detector. In this regime 
the system may develop a bifurcation with two possible 
stable states; it may be manipulated to force it to the 
state, correlated with the state of a coupled qubit. In 
contrast to the Josephson bifurcation amplifier, one of 
these states has a zero amplitude, which simplifies the 
task of resolving the two states. Furthermore, the prop- 
erties of the detector are different from those of the JBA 
for similar parameters. In particular, the switching curve 
may be narrower than that for a JBA, that may result 
in a higher fidelity of the qubit readout. 
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